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Abstract 

We get the general expression for the stabihty criterion with respect to the internal kink under helicoidal 
discharges, if one assumes that the helicity is finite and the aspect ratio is large {hx -C 1). 

We show that the helicity relative to the cylindric geometry has a stabilizing effect. This effect can be 
greater than the toric effect, when the helicoidal pitch is smaller that the toric radius. 

Resume 

On etablit I'expression generale du critere de stabilite relatif au kink interne pour les decharges helicoi'dales, 
dans I'hypothese d'une helicite finie et d'un grand rapport d'aspect {hx <C 1). 

On montre I'cffet stabilisant de I'helicite par reference a la geometrie cylindrique. Get effet peut etre plus 
important que I'effct torique, lorsque le pas helicoi'dal est plus petit que le rayon torique. 

Introduction 

La stabilite du kink interne a deja ete etudiee (dans le cadre de la MHD ideale) pour differentes config- 
urations, en considerant comme un petit parametre I'ecart a la geometrie droite infinie. On a ainsi traite la 
geometrie droite finie et circulaire [], la geometrie droite infinie non circulaire [], la geometrie torique finie 
[]. Dans toutes ces situations, il s'agit d'equilibres possedant une invariance par deplacement. Dans ce cas il 
restait a traiter la geometrie helicoidale de pas fini (et grand), ce qui est I'objet de ce travail. Pour I'essentiel, 
la methode utilisee reste la meme que precedemment []; elle sera rappelee au fur et a mesure de son emploi. 
Precisons que nous nous contenterons d'etablir I'expression de I'energie propre du mode, le traitement de la 
couche et le calcul du taux de croissance restant identique a celui de la configuration torique (verifier ^^p). 

Les coordonnees helicoi'dales 

(Indiquer les notations de Mercier.) A partir du systeme direct de coordonnees polaires r, ip, z on forme 
le systeme : 

r = r u = hip + £z v = —lip + hz 

ou les constantes h et £ sont normalisees par h"^ + 1"^ = 1 (avec > 0; ^ > ou ^ < 0). Chaque helice [H) 
definie par r = cte, v = cte sera un lieu d'invariance pour I'equilibre {z est normalise par tq). Sa transformee 
rotationnelle est j. Elle a pour periode en z 27r|, en u {H) est droite pour £ positif, gauche pour I 
negatif (axe magnetique). La necessite d'introduire des precisions (c/. Annexe A) provient essentiellement 
du fait que le systeme r, u, v n'est pas orthogonal. On utilisera une base de vecteurs orthogonaux, obtenue 
en remplagant Vu par un vecteur U = \7u — K\7u tangent a H . ha difficulte est alors en partie reportee et 
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se retrouve dans le fait que U n'est pas un gradient : on a rot [7 = ^1^^^:^ . Pour simplifier au maximum les 
calculs, il y a interet a efFectuer deux autres modifications : 

1) renormaliser la variable r pour simplifier la metrique. On posera x = |Vt;| dr, ce qui donne : 

2) introduire des coordonnees et des vecteurs complexes : 

^ X ~ iv X + iv -± -± 

S = — — t — — S = \7x + iVv 



T = Vx — i Vu 



Equilibres helicoidaux 

Dans le systeme de coordonnees x, m, w, I'equation du flux poloidal s'ecrit : 

df 1 dp 
dF^lPdF 



(1) v'LF + 2heU'f{F) + f% + l,%^0 



on 

Solution Generale. Le champH se mettant sous la forme : 

(3) B = f{F) U + U xVF 

on choisit I'unite de longueur pour que I'axe magnetique corresponde a r = 1 (x = 0) et u = 0. 

On supposera dorenavant que \x\ et \v\ restent petits dans la zone concernee |a;| ~ \v\ = 0(e) 1. Ainsi 
U, V et r resteront voisins de I'unite. L'ordering est standard. On a : a; et w d'ordre e, f d'ordre 1, F d'ordre 
e^ (done VF et Bp d'ordre e) ^ d'ordre 1, ^ d'ordre 1. Les derivations en a; ou u abaissent les quantites 
d'un ordre. Au point de vue signe, / est positif, F est decroissant a partir de I'axe. Les lignes magnetiques 
sont des heHces gauches pour g > (c/. Annexe D). Une perturbation a fc > est une helice gauche. 
Finalement ^ est negatif lorsque les lignes de champ tournent dans le meme sens que I'axe magnetique. 

On pent determiner aux deux premiers ordres (Annexe B) des families d'equilibres a "cercles" decentres 
(ce ne sont pas reellement des cercles mais des ellipses parce que x et w ne sont pas exactement des longueurs). 

\x\h<l 

[(Autre ordering, cf. Annexe B.) En fait j^j^ ^ ^ Mx suffit (c/. les termes correctifs de I'equilibre) d'oii 

ordering des termes non circulaire de I'equilibre.] 

A partir de I'expression du decentrement, on pent deduire celle de la modulation du champ poloi'dal le 
long d'une surface magnetique. Aux deux premiers ordres on a : 

\Bp\ ^ \Bp\o-[l + \h^pcose] 



■^Le champ de reference B = U comporte un courant (sans force). On pourrait prendre un champ de reference sans courant 
(par exemple B = Vu) mais la presence d'une composante poloidale impliquerait les choses. 
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oup et 9 sont les coordonnees polaires de la section et 



X + l= [ z^dz{Q^ - {4:5Q}+ I 2zdzP. 
Jo Jo 



On a. q = z = pi,qi et pi etant les valeurs de p, q et p sur la surface consideree = P, Q = 

Mqi =5. " 

Pour la suite, on introduira le coefficient de modulation "reduit" defini par : 

\ = \ + - + 5\\~X= j z^dziQ"^ -l-i5{Q = l)}+ [ 2zdzP 
4 Jo Jo 

A est nul si p et g sont constants. 

Coordonnees intrinseques 

On peut mettre le champ d'un equilibre donne helicoidal sous la forme 

B = V[ip- q{F) x] X VF 

ou X est independant dc (p ct multivoquc (0, 27r) [Hamada]. Dc U = Vw, — KV on tire B = Vu x VF + fU — 
KV X VF. Done on doit avoir [si u = ip] — V(<7x) x VF = fU — KV x VF. En multipliant par U : 

BpW\qx\=fU^ + KV-Bp = J = BW<p. 

Dans le Jacobian J, KV ■ Bp est du 2'^ ordre devant fU"^. [Le chouia fU^ correspond au fait que u n'est 
pas au depart exactement egal a la coordonnee d'Hamada (p telle que x = ~ Q^-] Dans le systeme de 
coordonnees intrinseques F, x on a : 



\dip q dxj 



on a dr = y^'^vy^Vx ^^^'^ VF • Vt^c x Vx = VF • [/ x Vx- D'oii dr = --^dF- djidx- Si on remplace F par 
H defini par — j dF = jidp,, on obtient : dr = y ^d/idud. L'element de longueur de I'helice est done 
I'element d'aire orthogonale aU est ^ ^idji dx- 

Dans le cas des equilibres a cercles decentrcs, on peut obtenir aux deux premiers ordres la formule de 
changement de coordonnees x,u,v Mj'^iX- H est dc la forme : x — iv = A(/i) + p,e~^" [w = a;(/x, x)] 



da = ^ = ^fid^dx, on a : 



on en deduit dxdv = (1 + A' cos ui) iJ,^dijdx- [H faut passer de s,t a z,z, z = ne^x] Comme on a 



(1 + A' cosw) ^ = ^T/2^w = x+ (/I'A + 2/i2 - nsmx 
2S = x-iv = n e-^^ + ^ {h^\ + 2h^ - £2)(e-2'x _ i) + a(m) . 



Rappel du probleme 
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Comme pour Ics configurations prcccdcmmcnt traitees, un petit paramctrc e (Ic rapport d'aspcct inverse) 
caracterise I'ecart de la configuration consideree au cylindre circulaire infini. Pour e nul, I'instabilite est 
marginals, et il s'agit de determiner, a I'ordre significatif en e, I'expression du reservoir d'energie 5w. 

Comme dans les autres cas, cet ordre est en par rapport a I'energie - positive - associee a ime pertur- 
bation arbitraire. Puisque nous ne cherchons pas a obtenir la perturbation elle-meme, nous pouvons profiter 
des proprietes variationnelles de 5w pour limiter la recherche de la fonction d'essai a un developpement en 
e moins pousse. Si la perturbation minimisante exacte est ^, si la fonction d'essai est de la forme ^ + i5^, 
I'erreur sur 5w sera 5w{5^). Cette erreur doit etre d'ordre superieur a e^. Cette condition sera remplie en 
imposant a I'erreur ^ deux contraintes : 

1) d'une part les composantes non contraintes de I'erreur 5£, seront d'ordre au moins egal a e^^o (soit 
au carre) 

2) d'autre part les composantes de I'erreur qui entrent dans 5w par une contribution en seront d'ordre 
au moins egal a e^. 

En d'autres termes, si ^ satisfait I'equation d'Euler F(^) = 0, on imposera a la fonction d'essai ^ + ^ de 
satisfaire + C) = O(e^). 

Une methode de resolution dite "alternee" est exposee dans 1' Annexe C. Une methode plus specifique 
dite "directe" est exposee ci-apres. 
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La methode directe 



Son application est limitee au kink interne. Elle consiste a construire une fonction d'essai pour le principe 
d'energie qui soit d'emblee exacte jusqu'a I'ordre requis. Ccici c'st possible parce que jusqu'a cet ordre elle 
est essentiellement independante de I'equilibre et qu'on peut done s'appuyer sur le cas simple oii le champ 
se reduit k B = U. 

Cette fonction ^ sera obtenue commc la somme de trois parties : ^ = Co + + ^2 (I'indice rappelant le 

mode dominant dans chacunc dc ccs parties). 

1) La partie principale ^ (qui sera determinee la premiere) s'annule au-dela de la surface resonante ; en 
dega, elle est strictement independante de I'equilibre (ainsi d'ailleurs que ^o)- 

2) Ci coincide a I'ordre principal avec le kink interne du cylindre infini circulaire. 

3) Co est une contribution simple de I'harmonique m = (une composante purement azimutale et uniforme). 
[Co est inutile puisquc C est defini a un multiple de B pres. Ah, mais si, B va changer. Mais alors, il faut 

I'introduirc quand on change B.] 

4) C2 ne contient que I'harmonique m = 2 (mais Ci en compte aussi). II est plus petit d'un ordre que 
I'harmonique m = 1 de Ci • C2 est adapte pour que C satisfasse les conditions de continuite sur la resonance 
et les conditions aux limites. II depend entierement de I'equilibre particulier. 

5) Chaque partie satisfait approximativement I'equation d'Euler, avec un residu qui est en valeur absolue 
d'ordre (par rapport a B^ aussi bien dans la direction azimutale que meridienne. 

6) Ci et Co seront exprimes en coordonnees helicoidales, C2 en coordonnees intrinseques. 

Le tenseur de deformation 

On utilise I'identite algebrique generale suivante : 

Q + J X C = -B • £> - V(C • B) 

oii le tenseur D est defini comme suit : 

D = VC + V^C--fdivC 

(en coordonnees cartesiennes : Dij = ^ + ^ — 6ij J|^). La construction de C consistera essentiellement 
a reduire I'ordre de D (le terme de pression ne presentant pas de difHculte). Ainsi </> sera essentiellement 

egal a — C • B. 

On appellera tenseur meridient Dm I'ensemble des composantes purement meridiennes de D et tenseur 
azimutal I'ensemble des autres composantes {D = £)„ + £>„). II sera necessaire de choisir C pour abaisser 

Du a I'ordre ct D„i sculement a I'ordre (il sera plus simple de dire qu'on annule D^ exactement). Pour 
cela nous calculerons d'abord les composantes de D sur la base S,U,T k I'aide de la formule generale {A et 
B quelconques) : 

A-D-B ^ ^V(C ■ B) + BV(C • A) - div [C ^ ■ + C • v4 X rot S + C • S X rot A . 
Nous poserons : Cs = C • ■S, Ct = C ■ C« = i C • 

AVC • B + SVC • ^ - ^ • B div C . 
Appliquons la formule de la trace (valable pour une base orthonormee) : 

^Ai + divC = 0. 
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[A reecrire en coordonnees normees et en base normee.] EUe donne ici 

U-D-U S-D-T 



(car XDX + VDV = SDT). On obtient pour L»„ : 

2 tt2 , 



+ div^ = 



UDU = 2k^U^^u-div{^U 



-SDU 
k 



-TDU 
k 



dS 



2M 



2M . 



Pour Dm on obtient [rot 5 = 0, rot T = 0, rot 11^=0]: 



SDS = 2V^{^-L^s 



2y2 = ST 



d^t 



dt 



TDT = 2V^ 
-SDT = +h^UiCs + &) + 2ev''ia- 
En ecrivant ^ = S + T — ikU ^ on explicite U-D-U sous la forme : 



U -D-U = k'^U'^£,u-U 



\- 




d 




dS 




^dt 





[/2 



div 



\2V^ J 



Faux : 



div^ = 



— +iK— 
dS du 



2^2 



\2V^ J 



d d' 
- — iK— 
dt du 



2^2 + « ^» • 



La fonction d'essai (^i et ^o) 

^1 sera clioisi parmi les solutions exactes (a trois ordres) de I'equation d'Euler relative au champ Bq = U. 
Dans ce cas (qui est sans force mais pas sans courant) la composante ^ ■ ?7 ne joue aucun role dans I'equation 
d'Euler. Nous pouvons done la fixer par la condition : 

^^7 + (/. = 

ce qui permet d'ecire I'equation d'Euler : 

U-D = 

(ce n'est pas une restriction reelle : la solution generale s'obtient a partir de la en prenant (p = - U, 
puis en redonnant k ^ - U une valeur arbitraire). A I'ordre le plus bas (et sans ecrire le facteur e'*^") nous 
devons imposer £^ - X = 1, ^ ■ V = —i, c'est-a-dire ^ • S* = 2, ^ • T = 0. D'autre part a cet ordre (j) (et 
par suite S,iU = — </>) est necessairement nul ; sinon Q* • G„ contiendrait une constante positive et 6w ne 
serait pas minimum ; cela correspondrait a la presence d'un mode m = de compression qui se propagerait 
necessairement jusqu'a la fronticrc. Par centre, un mode m = de rotation (sans composante normale) est 
acceptable jusqu'a la resonance : il sera decrit par avec une valeur de ^„ appropriee est une valeur nulle 
de (j) (cette decomposition en et est indispensable pour clarifier les calculs). 
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Revenant a ^i, il doit d'abord satisfaire I'equation d'Euler pour Bq, c'est-a-dire 



U-D = 0, 

ce qui donne [0(£^^)] : 

U-D-U = S-D-U = T-D-U = 0. 

Pour les elements meridiens de D, compte tenu du comportement de a I'ordre le plus bas, ils sont 
automatiquement d'ordre un. Mais il est necessaire qu'ils soient d'ordre deux pour que soit la solution 
generale a cet ordre. Puisque K est d'ordre un et nul a I'origine, on obtient pour SDS et TDT les conditions : 

On verra que cet ensemble de conditions fixe suffisamment. II reste a considerer SDT. D'apres son 
expression, il sera d'ordre deux si + k'^ = 0. Cela s'obtiendra en adjoignant a Ic tcrme forme 
purement d'une composante azimutale (autrement dit contient la modification de ^ ■ U qui etait a priori 
arbitraire avec Bq, et qui est appropriee pour un champ quelconque). [Soit div^ = sur I'axe.] Done pour 
1*0 : 

t ^ jku c ih U 



Le developpement limite de la solution 

La solution doit satisfaire les trois equations : 

UDU = , SDU = , TDU = 
et les conditions "initiales" sur I'axe magnctique : 

= = 5. = „. f = 0, f = „. 

Bien qu'elles ne determinent pas completement ^i, nous verrons que ces conditions suffisent a fixer les 

quantitcs qui joucnt ultcricurcmcnt. 

Une fagon simple de proceder pour etudier la solution consiste a eliminer ^5 et et a ecrire une equation 
pour ^„ de la forme : 

Les conditions initiales reportees sur ^„ donnent : 

-df-^ -ds-^V'^J' Is^-^- 

II est evident qu'il existe pour ^„ un developpement convenable (en commengant par ^5%), les termes en 
S^ et (pour p>3) pouvant etre choisis arbitrairement. 

Cependant, pour obtcnir les coefficients qui nous intcressent, il est un peu plus court de calculer directe- 
ment ^5. On verra qu'il est suffisant de calculer sur I'axe : 

i±SDS=^-2— 
2 dt dtdS dx 
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(compte tenu des proprietes de ^s, L = kK) 



-±SDT = +2h''^+h'U^+2k'V'^ 
dS ax db db 



j2 t 

L'inconnue fondamentale est done g^^- 



Posons 



Les trois equations liant S,u, X, Y s'ecrivent : 

(C) ^^^^ -2k'^ £ 

iUO) = 0, 1^(0) = 0, X{0) = 2, (L(0) = 0)). 

Unc elimination presque complete de ^„ entre (A) et (B) resulte de la quasi-permutabilite des operateurs 
Ds et Dt : 

On multiplie (A) par + L, (B) par ^ — L et on soustrait : 

On constate maintenant que (C) et (D) contiennent les memes derivees : ^ et qu'on peut separer 
par combinaison et ecrire facilement en ^5 et : 

^=^+^'[Wd^V^+^+2Ud^)'^'md^^^-[TPd^-^U 
rp^ n-d^t,^[V^ d U I dg\ I df ( dL P 

^-dS+^'[Wdiv^-^+md^)-^'md^ + ^-[-TPd^-^u 

et ^ sont les bonnes inconnues puisqu'on fixe pour t = et pour 5 = 0]. Nous avons besoin de 
deriver (E) par rapport a S. En fait, il sufRt de le faire sur I'axe en substituant a ^s, ^t, ^u, ^ 
les valeurs deja connues. Pour ^ et on les tire de (A) et (F) : 

^ = 2/ ^ ^ ^ 
dS dS dx 



On obtient (en faisant U = V = 1 quand ils ne sont pas derives) : 
= 



1 (Pis] , d fV^ d U 1 df\ 1 df dg f dL 



dx[^U dx ^ 2U dx) 4 dx dx [^dx ^ ^ 



2 dSdt 

[c'est bien ga la methode la plus rapide (pour abreger, poser Ds = — L, Dt = + L)]. 

4 dt dx \U dx 2U dx ) 4 dx dx 

d fl-V^ dU ^ d 1 ,1-1/2 df\ ^ 2 , dL 
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',^SDS-h'4- 

4 dt dx 



UI--5U 



+ t — 
dx 



[U'^ -l)[-U -5U'^)+2U 



?1 (--d-'u^ 

dx 



k'^{l-5'^) + k^{2.-5) 



1 d 



l-i, = l-r^C/^ = ^(f/^-l) 



= - — SDS + h^(--25 + 5'^ --\+k'^(-25'^ - — + — + ^8-1 



6^ 



On a ^5 = - L) C„ avec 



2k 



2U_ 



dS 



-L 



Ug \dt 



+ L 



dt 



+ L 



Yl 
Ug \dS 



L 



II sufEt de remplacer ^„ par {-jg — yi'^s 



2 dSdt 



d U 
di 



d^ u d r 1 df 

d^V^^di (2Ud^ 



Par substitution des expressions U, V, f, g, L on obtient : 

1 d^^s ,.(-. . . 9 



2 dSdt 



2(5 + J' 



4 dt 2 dtdS dx \ 2J ^ ' V 4 



d_ 

dS 



, 2M 

SDT = - 4 -5) + h\3 - 26) . 



Si on inclut ^0 dans la perturbation, on doit ajouter sa contribution a ^ SDT. 
De = ^ ^ e**^" on tire : 

5£>o T = - 2 A;2 y2 ^„ = + 2 V2 e'^"" 



4,SDoT=-4e^h^ = -kH\ 
dS 



Au total de ^0 et ^1 



dS 



SDT = -fc2(2 - Sf + /i^(3 - 25) 



On obtient la meme valeur pour ^ div^ = +^ du fait que le SDT total est d'ordre deux 
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Valeur residuelle de 

Puisqu'on annule les elements de D qui ne sont pas purement meridiens, seuls peuvent intervenir dans la 
valeur residuelle de BD les elements Dss, Dtt et Dst- Nous allons preciser leur contribution. Etant donnc 
qu'ils ont ete reduits a I'ordre final, et qu'ils seront multiplies par Q* avant d'etre integres dans SW, ils 
peuvent etre decomposes harmoniquement sur Tangle polaire de la section, en ne retenant que I'harmonique 
homologue a I'harmonique dominante de ^. 

Si nous definissons, a I'aide du vecteur unite normal e„ et du vecteur unite tangent a la section et les 
vecteurs : 

M = Bn + iet et N = — iet 
et si on appelle a I'angle {ex, e„) on a : 

M = Se-'"' N = f e'" 

d'oii Ton deduit aisement : 

4:ietDen = e"^'" SDS - e^'" TDT 

-4 etDet = e'^'" SDS + e^'" TDT - 2 TDS . 

Les termes significatifs des residus de SDS, TDT, SDT sont du 1'"' degre en 5 et i et varient done angu- 
lairement comme e~"* et e*". D'autre part les termes qui contribuent au mode m = 1 dans Bm Dm sont les 
termes en e~'". On voit done que TDT ne contribue pas, que SDS contribue par son terme en t et TDS 
par son terme en S. 



La contribution de 

On a note plus haut que pour un champ Bq = U, la composante ^ • f7 ne contribue pas a G. Autrement 

dit 

UDi^o) = V(eo • U) 

[^0 est-il parallele k U ou k B 1 (Ici a U, mais ne pourrait-on le prendre parallele a B ?)] Par un champ 
quelconque, on aura : 

G(Co) = BD(eo) - V(6 • (-(i3V)-MivCoVp) 
= fU-D + BmD-^{f^o-U) 
= {f-fo)U-D + W[{fo-f)^o-U] + Bm-D 

on fo est une constante arbitraire. 

Par la reciprocite, on a B^ ■ D - U = B^ V(^o • U), done 

Bm-D=^Bm V(^o ■U) + Bm-Dm. 

On a vu que dans D^, Co modifie seulement SDT, en le faisant passer a I'ordre deux et en changeant un 

peu les coefficients. 

Les deux termes (/ — fo)UD et Bm V(Co • U), sont negligeables : ils sont en effet d'ordre f £, dans 
la direction meridienne et £^ / ^ dans la direction azimutale (a comparer a £^ / ^ et £ / ^ comme on le montre 
plus loin). On retiendra done seulement : 

G{^o) = V [{fo -f)^o-U]+ Bm Dm{^o) - (BV)-i div Vp . 
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Le calcul de SW 



Soit I'expression : 

SW=^JdTQ*-G 

que nous appliquons a la fonction d'essai + dans le domaine limite par la surface resonante {qk = 1) 
on a 

G = BD + V[(/o -f)^o-U]- V(a • B) - (SV)-i div^ Vp . 
A I'ordre deux, B ■ D se reduit a Bm Dm{^o + ^i)- Nous partageons I'integrant en deux parties : 

G = Gi + G2 

avec 

Gi = BmD^-{p-pi)Vv BVr? + dive = 
G2 = -V(6 • B) + V[r?(p - pi)] - V[(/ - /i) ^0 • U] 

Pi : valeur de la pression a la limite. 

[On a remplace /o par /i, valeur de / sur la surface resonante.] Nous allons montrer que les termes de 
Gi sont d'ordre deux. G2 ne comporte que des gradients exacts qui donnent une integrale de surface par 
integration par parties. Indiquons d'abord que I'ordre significatif pour la composante meridienne de BD, 
{B-D)m, est I'ordre deux, (e^ /^) qui se combine avec I'ordre zero (/^) de Q^. Pour la composante azimutale 
de BD, soit {BD)u, c'est I'ordre un (e/^) qui est significatif, parce que est d'ordre e/^ (comme on le 
voit dans la methode alternee). 

Done les tcrmcs d'ordre superieur a ccux-la sont negligeables (comme on I'a dit pour ^g)- 

Dans Gi, B^ est bien d'ordre e"^ f Quant au second terme, il faut noter que div^ est de I'ordre de 
^ et en outre constant a I'ordre principal, done i] (en |) est aussi constant a I'ordre principal et Vr/ reste de 

I'ordre de j ; comme p est en e^/^, on a bien le meme ordre final. 

On integre par parties le terme contenant G2 et on obtient : 

SW^^J dTQ*-Gi + ^ J dSC -nBViii- B]. 

On va calculer separement le terme de volume et le terme de surface. 

Evaluation des termes de volume 

Soit 

SW=^J dTQ*^D^B^^^- j dTQ*^-Wv{p-Pi). 
A I'ordre principal ou ^j^ est constant, on a : 

avec 

a ^ (ep + I eg) e:'-^'^- = 1^.1 S er^'^^ 

[Bm = S,T, = 'S', r, Qm — S,T, Qm B„i D„i —], cc qui donne (par unite de longueur), pour le premier 
terme: 



5Wi 
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En substituant \Bm\ = ^ et {Dm ■ eg e'^) = — fc^ [i a + (a + 2^) eg] on obtient, apres une integration par 
parties sur g : 

( fp^ / 1 \ fp^ /I \ / 1 

|2 J-2 7^4 I „, / „3 j„ ( ^ 1 \ I fl / „3 .■ ' . \ / i 



Pour le terme de pression, on a : 

-6W2 = drQ* -Vvip-Pi) 



dr ip - pi) div {C Vr]B- CB Vrj) 
dr {p — pi) div (^* div ^) 



Or on a etabli 



= ljdT{p-pi) [div C div ^ + C V div . 

S-D-T U-D-U ^2 ito 
'iivC = u-2 +^ 

et puisque div^„ ~ div^o = -/i^e'*^", div^^ ~ 2/i^e''^". Enfin 

CVdiv^^^j^^vf ^^^2 cT^j =4fc'/3- 



Au total, on a : 



-SW' = er[^+^)lJdr{p-p,)4k\ 



Calcul du terme de surface 



Ce terme s'ecrit 



6Ws 



dS^^-nB\7{^i.B). 



Cette integrate doit etre effectuee sur la surface magnetique exacte oil le mode m = 1 est resonant {kq = 1). 

L'operatcur i?V fait croitre d'un ordrc Ic second factcur, puisque rharmonique m = 1 est dominante 
dans ^1. [D'accord mais avec quelle fonction-test ?] Et de ce fait, I'integrale prend automatiquement une 
valeur d'ordre deux. [Laval : V et ^ ou Pour tirer parti correctement de cette propriete, il suffit 
d'efFectuer cette partie du calcul en coordonnces intrinscques de Tcquilibre. Plus prcciscment, on decompose 
harmoniquement en x chacun des deux facteurs (a I'ordre un). Dans le systeme /x,w,x on a : 



BV = j{ik + 



ld_ 

q dx 



dS = 



dT\v^i\ /mIVmI 



dfi 



J 



dudx 



5Ws, = ^-j dudx (mCi* Vm) [ik + (a • B) . 

Par unite de longueur en u (avec kq — l) 



5Ws, = ^ I rfx[M^rVM] 



dx 
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B et ^- Vi^ sont Ics 2 quantitcs continues (a i?V pros) sur la resonance.] L'harmonique m = 1 ne contribue 
pas. A I'ordre considere, interviennent seulement to = et m = 2. L'harmonique to = ne joue pas (il n'y 
a pas de mode de compression par construction, et le mode de rotation a ete elimine par ^o)- 

II rcste done seulement a calculer I'apport du mode to = 2. On ecrira pour cela ■ V/U et i^i ■ B 
d'abord dans les coordonnees s, u, t, puis on effectuera le changement de coordonnees en fi, u, x- On pourra 
d'ailleurs faire d'emblee t = 0, puisque seuls les termes en 1, S, contribueront significativement a e~'^^^. 

On peut ecrire (c/. Annexe C) : 



u_ 
IP 



X = ji^- V/i = i 



dW 
dx 



i^-B = i^-Uf + i^-Bm = 



dxJ 



k dx 



(expressions valables aux deux premiers ordres ; W et n'ont pas de composante en to = 0. Par centre 
i^- B en possede une d'ordre deux, qu'il faut supposer absorbee dans ^o)- Pour determiner W{S) on utilise 



W{S) = 2S+{2f -h'^)S'^ . 



D' autre part : 



2S 



2 

1 - 
jj e 



^ 2ih 
k 

2h^S 

9 

-iX 1 



S + 2S^ 



k 



\h^\+2h'' -e) + 



3h^) 



ce qui donne par substitution : 



W 

X* 
+ i<j) 
SWs, 



= lie 



2h£ 



A + 1- 



2M 



lj.CWii = n'^h^e^'^ (^X+^^ pourm = 2 



-i^.B = -fkjjL^ h'^e 



-2ix 



3,, ,x hif^ 1 
2(A + 1)-^ A+- 



h£ 



2(A+l)-y(A+, 



1 



^4 



pour TO 



Complement de la solution 

La solution complete doit satisfairc les conditions aux limitcs, ainsi que les conditions de continuitc sur 
la surface resonante : continuite de Q ■ \/F et de B ■ G (pression totale perturbee), autrement dit de B ■ V0 
pour une solution de I'equation d'Euler. II est done necessaire d'ajouter a + Ci une troisieme partie ^2, qui 
dcpendra de rcquilibrc. Si on decompose la solution en liarmoniques de Tangle intrinscquc x- I'liarmonique 
dominante to = 1 n'a pas de condition de continuite a satisfaire et sera nuUe au-dela de la surface resonante 
en raison des conditions aux limites. Pour l'harmonique m = 0, Q ■ VF et de B ■ G sont nuls sur cette 
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surface, done aussi au-dela pour la meme raison que m = 1. Ainsi ^2 ne contient qu'une harmonique m = 2, 
de part at d'autre de la surface. Son amplitude est d'ordre s devant ^1. II suffit par suite de la calculer a 
Fordre principal meridien, suivant la methode de I'Annexe C. La contribution de ^2 a SW apparait comme 
un terme de surface complement aire. On a en effet (en n'ecrivant pas ^0) 

SW = SWi^i) + 2 6) + 5W{^2) 

avec 



puisque le terme de volume du couplage est d'ordre superieur (a condition d'ecrire I'integrant dans le sens 
convenable). On doit determiner pour ^2 la solution interieure et la solution exterieure ^e- Elles sont 
fixees par Ics conditions a rcmplir sur la resonance (ainsi que les conditions aux limites et les conditions de 
passage sur une eventuelle surface resonante m = 2). On doit avoir sur la resonance : 



+ 01 = 



{X = /i^V/i, (}> = {Xi et ctant les composantcs m = 2 dans ^1). Les conditions sur I'axe et sur 
la frontiere determinent les relations (pi = AXi et cpe = B X^ {A et B dependent du profil de courant). 
En calculant et en ajoutant les differents termes de surface, on montre aisement qu'on obtient le total en 
remplagant dans I'integrale de surface de ^1 seul le produit ^1 (pi par <Pi — ^i^e, c'est-a-dire : 



ce qui donne : 



avec 



AXi]\ 



BXi] 



B-A 



SWs = irkf 



B-A 



h-AXi] [^i-BXi] 



(A + 1) 



he 



x + 



= /x/i(A + 3/2) 



A et B sont calcules dans I'Annexe C. Pour m = 2 et fcg = 1, on montre que : 

"1 M 



A on B = kf 



Pour un courant plat on a : 



d ^„ hi 3 
— £VuH 



1 



d 



4 ^V/i djjL 



1 



4 ^V/i d/Lt 



3 
"4 



on posera done 



A = 
B = 



kf 
kf 



k -h' 
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En remplagant d'autre part A par A 



on obtient 



A + a I 



M 
k 



+ b 



1 + a-b 



Regroupement des resultats 

Pour condenser les formulcs nous mettons a part le facteur de normalisation de 6W qui est egal par unite 
de u a 7r/^fc**/Lt"^^^. Et nous posons 



hi 



Nous avons obtenu SW comme une somme de 3 termes : 
1) un terme de courant : 5Wi = aj+/3(|^+r) avec 



5 
T 



^0 



4 / z'^dzi —-1 
kq 



1 



1 



/3 = -{i-6r + 



fc2 V4 



2) un terme de pression : 6W2 = avec 



P 
Y 



- V 

JO 



zdz 



P-Pi 

■"mi 



3) un terme de surface 

6W3 



1 

1 + a + b k"^ 



A + a A 



3 
4 



avec X = T-6S + P. 

Get ensemble se met sous la forme : 



A-2 



'+2 



(l_2<5)-(l-<5)2- 



T 
2 



1 + a-b 



bX^-abiX+--d 



Si on fait h = 0, {S = 0), on retrouve le cylindre : 



sw = -{x + ^ 



P-Pl 
52 
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,4 

Si on fait ^ = mais ^ ^ (helice infinie, Stellarator) on obtient : 



5W = -{l-5f 



T + P+[ -^-5\S 



(T + P) 1 
stable pour S > ^ h - 



efFet stabilisant pour 5 > 0. 

Si on fait £ = 0, on retrouve le tore : (5 = 0) 



SW 



T + P- 



S 



[1 - ] + + 



T ^a-b{T + Pf ~ab{T + P+iy 



1-a-b 



Si on fait 6=1 (equilibre sans courant) 



Discussion — Interpretation 



La formule contient deux sortes de parametres : 

1) des parametres caracteristiques du profil de courant (et de pression). Ce sont les meme qu'en geometrie 
torique, a savoir : S, T, P, a, b 

2) deux parametres de forme, un pour I'helicite |, I'autre pour la periode normalisee de la perturbation, 
(5 = M = d S)^^, on periode des lignes magnetiques resonantes, rapportee a la periode hk des lignes 
magnetiques d'un champ sans courant (lorsque le courant axial total du au champ poloi'dal compcnse le 
courant dn k B = U). Ce champ correspond a 5 = 1 et constitue en fait le veritable champ de reference (le 
champ en I'absence de plasma). 

D'ailleurs on constate que pour 6 = 1 (equilibre sans courant), SW est (a un facteur pres) independant 
de I'helicite. 

1) Pour discutcr Ic rcsultat obtcmi, nous obscrvcrons d'abord I'cfFct dcstabilisant de la pression (a travers le 
terme b essentiellement) et nous examincrons cnsuite ce qui se passe pour P = 0. 

2) Nous supposerons le shear assez faible pour que a et 6 soient petits, ainsi que T et S. Nous retiendrons 
les termes du premier ordre pour obtenir : 

SW ^,3 3/^4 / 3 , 29 

oil on a fait T ^ 25 et a ^ 45 (shear faible). 

On verifie que le coefficient de p- est toujours positif. On pose done p- = (5^ ^ et on a la condition de 

{1-Sf 



stabilite : 



^ > 



m 



P3S^6^-IS + f,} 

d'oti Ton deduit que la condition est toujours remplie pour e > 1 (helicite des lignes magnetiques inversee 

par rapport h U ct courant principal inverse). 

Dans les autres cas, on peut s'assurer que la fonction F{S) decroit de +00 a quand S varie de zero a 
un, et decroit de +00 a zero quand S va de zero a —00 (helicite de meme signe). II en resulte qu'une valeur 

donnee de ^j- stabilise tous les S assez grands (en module) c'est-a-dire les modes de periode assez grande. 

.2 

Et la stabilisation est d'autant meilleure que ^3- est eleve. 
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Si on compare maintcnant avcc Ic cas de la geometric torique (quasi-marginale), on remarque que la 
stabilisation provient ici de nouveaux effets, a savoir essentiellement : 

1) du fait que la courbure (^) est remplacee par qui peut etre grand en pratique, meme pour une helicite 
reduitc parce que le pas de I'hclice peut etre notablement plus court que le pas de la perturbation resonante 
(des lignes magnetiques), tandis qua dans un tore ils sont necessairement identiques. [Torsion?] 

2) du role de 6 dans le formule, en particulier la stabilisation pour 6 > 1. 
Par exemple, pour 6 — )• — oo on a la condition : 

5 h? 1 + a + h 

[Si on pose Q = 45^ - (25 + T + P) et 5 = 1 + r? ga fait : 



5W 



Q 



1] 



l + 2r] 



T 
2 



(l+a + b)-' + +b{Q + 2Sf +ab(Q- 



2S- 



Du 2** degre en Q ou en ^.] 

Lorsque S > 1 + ^"'"4^^'^ le terme independant de ^ est positif. Alors on est stable si ^ n'est pas trop 
grand, ou si son facteur est positif aussi, c'est-a-dire essentiellement bX^ pas trop grand (on peut ecrire la 
condition |6A^| <). 

Annexe A — Coordonnees helicoidales 



Posons 



Vu ■ Vv Mir^ - 1) 



Soit 



X = yx U = Vu- KVv = 



= K{r) 



+ hr e 



dK 



M U" 



V = Vw = hcz eu 



Soit f/ = |f7| et F = |y| = ^ = \X\, X, U, V sont orthogonaux (mais non normes). {U 
rUV = 1). [Commutateurs : {U, V) = 0, {U, f ) = VV^ x V, (f , V) = rot {V^ U).] 
On a divC/ = 0, diY {U X) = 0, divF = 0, rotX = 0, rotF = 0, rotU = 2MU^U. 
On utilisera : !^ = ^ = -fU^V, K = M U^{r^ - 1), ^ 



S=^, t 



x+iV 



2hi^, VS = Vi = f . Avec 



on a : 



S + t V = i{S-t) 



d d . d d d . d 
dS^di'^^dV Jt^di~^dV' 

Soit S = X + iV , U, T = X — iV la, base complexe. Toutes les quantites suivantes sont nuUes : 

5^ f2, S-U, f-U, mtS, rotf , div(C/5) , div{Uf) , divC/. 

On a : S-f=2V^, plus \U\^ = , lotU = 2h£ C/^ u. [SVS = 0, TVT = 0.] On a encore Sxll 
etfxU = iUf 



-i US 



5 V/ = 

TV/ = 



dt du 



SxT = 2i — U. 



Pour les quantites de rcquilibre = 0, et pour les perturbations (^ = e''^" ^(r, V)), ^ = ikf et UVf = 
i k W^f. (Donner le tableau AX/B pour la base R, U, V.) 
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Annexe B — Equilibres a "cercles decentres" 



Comme en geometric torique, les equilibres a fronticre circulaire ont des surfaces magnetiques dont les 
sections meridiennes sont des cercles decentres, aux deux premiers ordres (du moins formellement, dans les 
coordonnees x ct V qui ne sont pas de vraies distances). [C'est vrai aussi pour un cercle reel, car I'ellipticite 
est petite.] A I'ordre principal, le flux F = Fq satisfait : 

Posons X = pcos9, V = psinS. 

A cet ordre, les solutions a frontiere circulaire ne dependent que de p. Soit Fo{p) I'une d'elles. On a : 



(1) 



^ + l^ + 2.W) + /^(i^o) + ^(Fo) = 0. 



dp'^ p dp 



dF 



dF 



A I'ordre suivant, comme tous les termes correctifs qui constituent le second membre sont en cos^, il en est 
de memo pour la modification a apporter k Fq. II est commode de chercher I'ensemble du flux modifle sous 
la forme : 

F = F{p, 6) =Fo[p- A{p) cose] = Fo[a{p, 9)] . 

Le decentrement A{p), d'ordre ep, du "cercle" a = cte (en realite une ellipse) sera regi par une equation 
qu'on obticnt en ccrivant que F satisfait I'equation d'equilibre aux deux premiers ordres. Par derivation de 
I'expression de F, on a d'abord : 



1 d dF 1 d^F 
p dp^ dp p^ dO"^ 



da^ 



1 dFo 
a da 



■ COS0 



D'autre part, en prenant a pour argument de Fq dans (1), on a 

Nous additionnons ces deux egalites et nous deduisons de leur somme I'equation de I'equilibre exact (qui 
doit etre valable pour F a deux ordres). Nous obtenons : 



cos^ 



F^A" + 2F^A' 



dF 



1 



dp 



En retenant I'ordre principal (termes en cos^), en multipliant par pF^ et en integrant, il vient {Lq F 
F" + \F'): 



p^F['A[+epiF, 



Jo 



dp 



pF'^ + 4h£p'^fF' -2p- 



dp 
dp 



= 0. 



Nous posons F' = -p^ {F' est negatif) et ^ = Nous obtenons : Mqi = 5, = Q, ^ = P' 



Pi 



La modulation du champ poloi'dal sur chaque surface magnetique est liee directement au decentrement. Bp 

etant de la forme U x VF((t) avec a = p ~ A cos 6*, la modulation relative de \Bp\ est la somme de celle de 
\U\ et de celle de |VF| (a a constant), c'est-a-dire de Vcr, ou de Vp(l — A'cos0) (car la contribution de V9 
n'intervient qu'au second ordre). Or |Vp| est module comme |Va;| et |VF| c'est-a-dire comme V (I'ellipticite 
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dcs "ccrclcs"!), dc sortc qu'au total \Bp\ est module commc UV{1 — A' cos 61) ^ 1 — (p + A')cos^. Nous 
ecrivons cette expression sous la forme 1 + h'^Xp cos 9 en posant A' + p{l + Xh^) = 0. Soit 



X + l= [ 7Jdz{Q^-4SQ}- [ 2z'^dz[P']. 
Jo Jo 



Cette quantite interviendra dans la deuxieme partie. Les equilibres obtenus sont valables lorsque la pertur- 
bation de F reste petite devant Fq on constate que la modulation du champ poloi'dal est en h^. 

On aboutit a la meme conclusion par la minimisation de 6W. En effet la contribution azimutale positive 
Q* • Gu doit etre reduite a I'ordre etant donne I'ordre du couplage de Q* et G„ avec QJ^ et Gm- 

En geometrie helicoidale, (2) s'ecrit : 

U-G = BmV{^ ■ U) - /div [U U"^] + 2MU^£,- WF = 0{s'^) . 

Puisqu'il y a un arbitraire dans ^, on pourra supposer que ^ ■ U est d'ordre e^m, ou alors d'ordre mais 
constant a cat ordrc. Des lors tous les termes autres que la divergence sont d'ordre e^. (2) s'ecrit done 
div [$m U"^] = 0(£^) ou encore 

(3) ^ = X ^ + 0{s^) • U fait nul) . 

Pour faire le pas suivant, on reporte I'expression (3) dans la composante meridienne de (1), qu'on va pouvoir 
resoudre aux deux premiers ordres (I'ordre zero est I'analogue de I'ordre "cylindrique" usuel, I'ordre un est 
celui des couplages de modes et aussi I'ordre ou joue la pression) . 

Cette resolution va donner ^ aux deux premiers ordres et en meme temps </> aux deux premiers ordres. 
Bien entendu, a I'ordre zero, on peut choisir a volonte le nombre d'ondes m du mode qu'on calcule. 

On peut par exemple eliminer <p en ecrivant : 

(4) U-TotG = 0{e^) avec (3). 
En negligeant des quantites d'ordre on trouve : 



(5) -[/•rotG = div 



C/^V, ^^^^^ 



+ VW ■ rot Ju + Vp- rot [[/(BV)"^ {VW ■ rot)] . 



V J. 

II est aussi necessaire d'evaluer les quantites qui interviennent dans les conditions de continuity. Ce sont : 
Bm-G ^ BmV<P ^ fBVWX B^VW J_U_ 

JJ2 U2 ^ y U2 J+ JJ2 JJ2 ■ 

La composante de qui est egale a 2h£ ne joue pas a I'ordre principal dans I'equation d'Euler de W 
parce que c'est une constante mais elle joue dans I'expression de ^. 

A I'ordre principal, apres multiplication par ^ , on obtient : 



J_#_/J__d__ \ dW_ ^/hi l\.dW 
kf dx \kq dx J ^ dp \k qkj * dx 



Pour I'harmonique m, en introduisant 4 • V/i = — W 



, , M f / 1 1 \ M d ,^ ^ ^ 2M 1 1 
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Une fois resoluc I'cquation dc W, on rcvicnt a I'cquation azimutalc cn y introduisant la valeur dc (/) qui 
est maintenant comme aux deux premiers ordres. On doit done satisfaire I'equation de U ■ G avec un 
second membre, en tenant compte cette fois de tous les termes d'ordre e^, avec des valeurs approchees de ^, 
c'est-a-dire exprimccs cn W. Autrcmcnt dit la correction porte seulement sur div(^m [/^). 

Dans le cas du kink interne, il est necessaire d'aller jusqu'a ce stade (2 ordres meridiens, 3 ordres 
azimutaux) pour construire la fonction-test. Pour obtenir la vraie perturbation singuliere, il faut meme aller 
au 3^ ordre meridien (selon le meme principe iteratif). La formule (3) reprcscntc ncccssaircrncnt toutes 
les solutions aux deux premiers ordres. (C'est par exemple une consequence evidente de UDU = 0). En 
appliquant la formule Bm ■ D ■ U = on obtient (pour la solution particuliere) 

Au deux premiers ordres on peut remplacer i ^ VF par Bm . 



Annexe D — Sens de rotation des differentes helices 



1) Pour definir une helice on choisit un axe (oricntc) Oz, un rayon r et une helicite (avec un signe). [Sens 
direct ou retrograde ? (definir).] Sur un petit cercle de rayon r et d'axe Oz, on fixe une orientation (droite) 
pour Tangle ((/?), de telle fagon que le triedre er, eip, ez soit droit. On choisit un point arbitraire de I'helice. 
L'helicite fixe le rapport constant de I'accroissement de a celui de z le long de I'helice. Si elle est positive, 
I'helice est droite, si clle est negative, I'helice est gauche. 

2) L'helice {H) (I'axc magnetique) a pour helicite j. [Toute helice de I'equilibre.] Puisque on prend h 
toujours positif, [H) est droite pour C. > gauche pour (. < Q. Si r, y, z est droit, r, u, V est droit (le 
jacobien valant un). Sur [H) V = —£ip + hz est fixe. 

3) Le systeme intrinseque ii,u,x 6st droit {fidiidudx = y drdudV). Pour les lignes magnetiques, le rayon 
est /i, I'axe est u, Tangle est x- H ^ut done considerer pour reperer les helices, le systeme jj,, x, u qui est 
gauche. On a u — qx = cte sur une ligne. Done une ligne magnetique a q positif est gauche. [La periode 
d'unc ligne magnetique est Am = 27rg = ^ longueur 2tt ((Jj:) rapport d'aspect fc/j, -^.J 

4) Une perturbation rcsonante doit avoir Thclicitc d'unc ligne magnetique, done ctrc gauche si q est positif. 
Cela se traduit par le fait que A: a le meme signe que q. Son facteur de pcriodicite est d'aillcurs eJ'^'^'^^i^) , 

5) En conclusion, si ou Mq est positif, la perturbation a une helicite de signe oppose a celle de Taxe ; si 
^ est negatif, les helicites sont de meme signe. 

6) De meme pour les courants principaux U [2hi f U^] et U LF. Le second est egal k —U f 
II est de sens oppose au premier si Mq > 0. 
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